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1  Description  of  Research  Accomplishments 

1.1  Overview 

The  major  aims  of  this  grant  as  described  in  the  original  proposal  were  to 
1)  develop  improved  analytical  models  for  multidimensional  reactive  flow 
and  ignition  of  energetic  materials  2)  improve  computational  capability  for 
reactive  flow  and  ignition  problems  3)  transfer  new  algorithms  to  explosive 
design  codes  and  4)  document  advances  in  the  theory  and  computation  of 
reactive  compressible  flow  relevant  to  the  explosive  system  design  process. 

The  support  given  by  the  AFRL/Munitions  Directorate  over  the  past 
eight  years,  (from  this  grant  and  a  previous  three-year  contract)  has  enabled 
dramatic  advances  in  detonation  theory  and  numerics  for  application  to  the 
design  of  explosive  systems.  In  particular,  D.  S.  Stewart  in  collaboration 
with  J.  B.  Bdzil  of  Los  Alamos  National  Laboratory,  have  pioneered  the 
theory  of  Detonation  Shock  Dynamics  (DSD)  as  documented  in  numerous 
archival  publications.  The  DSD  method  uses  theoretical  analyses  of  the 
reaction  zone  structure  to  determine  how  the  lead  detonation  shock  moves. 
In  particular  one  determines  the  relation  between  the  normal  detonation 
shock  velocity,  the  shock  curvature  and  other  higher  derivatives  of  the  shock 
surface.  This  analysis  then  leads  to  simplified  methods  to  predict  the  shock 
motion.  The  methods  also  lead  to  ways  to  characterize  classes  of  explosives. 

The  University  of  Illinois  has  been  the  focal  point  and  originator  for 
many  of  the  most  recent  theoretical  and  numerical  advances.  In  particular, 
prototype  numerical  codes  have  been  written  for  DSD-characterization  of 
condensed  explosives  with  real  material  properties  and  for  DSD- wave  track¬ 
ing  and  burn  time  calculations.  DSD-technology  has  been  installed  in  major 
National  Laboratory  hydrocodes  and  that  activity  is  ongoing.  The  scientific 
predictions  of  the  DSD-theory  for  explosives  have  been  validated  by  direct 
numerical  simulation  (DNS).  New  engineering  strategies  have  been  devel¬ 
oped  for  explosive  design  based  on  this  relatively  new  science.  New  models 
based  on  continuum  mechanics  have  been  developed  for  modeling  the  ig¬ 
nition  of  energetic  material  from  its  realization  as  a  solid,  all  the  way  to 
a  reactive  gas.  Most  if  not  all  current  trends  in  national  defense  applica¬ 
tions  that  involve  insensitive  munitions  and  miniaturization  for  explosive 
systems  can  make  effective  use  of  this  technology  and  continuing  research 
developments  in  this  area. 
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1.2  Chronology  of  Activities 

In  what  follows  below,  we  give  a  chronology  of  research  activities  and  de¬ 
velopments  that  were  undertaken.  We  make  reference  to  the  papers  and 
reports  that  documented  this  work  and  in  Section  2.  one  can  find  abstracts 
or  summaries  of  those  reports,  with  more  details.  All  of  the  papers  and  re¬ 
ports  and  associated  computer  codes  have  been  delivered  to  Eglin  engineers 
throughout  the  course  of  the  grant.  Additional  copies  of  these  materials  are 
available  on  request. 

The  fall  of  94  was  effectively  the  end  of  our  previous  contract  with  the 
Munitions  Directorate.  The  status  of  DSD  -  theory  and  related  computa¬ 
tional  algorithms  had  significantly  advanced  by  the  work  on  that  contract. 
A  new  set  of  robust  numerical  methods  based  on  level-sets  had  been  iden¬ 
tified  for  propagating  DSD  shock  wave  in  three  dimensions  and  prototype 
codes  had  been  written.  The  effect  of  shock  acceleration  on  DSD-theory 
had  been  considered  for  the  first  time.  The  numerical  capability  to  work  on 
complex  multi-D  problems  with  embedded  boundaries  had  been  significantly 
advanced. 

However  most  if  not  all  the  developments  in  numerics  and  theory  had 
been  carried  out  for  an  explosive  described  by  an  ideal  EOS.  Hence  the  work 
of  this  grant  focused  on  theory  and  numerics  for  explosive  with  real  proper¬ 
ties,  described  by  non-ideal  EOS.  Also  work  was  carried  out  to  implement 
DSD-routines  and  transition  them  into  AF  and  Department  of  Department 
of  Energy  (DOE)  hydrocodes. 

1.2.1  Spring  1995:  Upgrading  the  Computational  Capability  of 
the  UIUC  Group 

In  the  first  few  months  of  this  grant,  in  early  1995,  a  significant  upgrade 
of  the  groups  computational  capability  was  achieved  through  the  purchase 
of  an  SGI  Power  Challenge  Server  with  1  Gbyte  of  RAM  (machine  name 
Messiah).  With  the  purchase  of  this  machine,  Stewart  no  longer  needed 
to  write  separate  proposals  to  the  NSF  supercomputing  centers  to  get  ad¬ 
equate  computational  cycles  to  carry  out  the  research  on  the  Directorates 
projects.  The  turnaround  time  on  computational  runs  in  the  group  dropped 
dramatically. 

Because  of  the  need  to  validate  and  to  test  DSD  algorithms,  the  hy¬ 
drocodes  (DOE)  MESA  and  (AF)  EPIC  were  ported  to  the  SGI  Power 
Challenge.  Both  codes  were  chosen  as  target  hydrocodes  for  the  installment 
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of  the  DSD  wave  tracker  (DKAPPA2D),  Changes  to  program  burn  algo¬ 
rithms  motivated  by  DSD-based  analysis  were  identified  through  collabora¬ 
tion  between  Stewart  and  John  Bdzil  of  Los  Alamos  National  Laboratory 
(LANL). 

1.2.2  Spring  1995:  Design  of  AXS  a  New  High  Order  Reactive 
Flow  Code 

In  order  to  carry  out  validation  studies  of  DSD  and  comparison  with  Di¬ 
rect  Numerical  Simulation  (DNS)  of  reactive  flow,  it  was  deemed  essential 
to  write  a  new  hydrocode  that  employed  high  order  methods  with  1990s 
algorithms.  Previous  attempts  to  measure  detonation  shock  dynamics  of  a 
DNS  in  converging  and  diverging  wedge  geometries  had  been  carried  out  by 
Tariq  Aslam  (a  UIUC  graduate  student  supported  by  this  grant),  Bdzil  and 
Stewart  with  the  use  of  (DOE/LANL)  CAVEAT  and  MESA.  Low  accuracy 
measurements  of  the  relation  between  the  normal  shock  velocity,  Dn  and 
the  shock  curvature,  n  were  had  been  made  with  MESA.  But  an  inherent 
difficulty  was  that  the  measurement  of  the  shock  curvature,  the  shock  ve¬ 
locity  and  the  shock  acceleration  required  construction  of  derivatives  of  the 
motion  of  the  shock  surface  and  lab  codes  like,  MESA,  CAVEAT  and  EPIC 
did  not  have  high  accuracy  algorithms. 

In  April  of  1995  the  new  code,  named  AXS  was  designed  in  a  collab¬ 
oration  between  Stewart  and  his  two  students,  Aslam  and  Xu.  AXS  used 
EN 0-schemes  that  were  4th  and  5th  order  spatially  accurate.  The  work 
was  documented  in  [6]  and  AXS  was  delivered  in  1998  to  Eglin  engineer 
Michael  Nixon.  Two  versions  were  written  corresponding  to  the  two  stu¬ 
dents.  Aslams  version  implemented  reactive  flow  for  an  ideal  equation  of 
state,  but  importantly  Xus  implementation  was  for  nonideal  EOS.  The  sec¬ 
ond  version  of  AXS  was  used  to  compute  deflagration  to  detonation  tran¬ 
sition  in  porous  and  the  nonideal  version  of  AXS  implemented  a  Glaister 
Riemann  solver.  This  work  was  documented  in  [3].  The  details  of  the  algo¬ 
rithms  can  also  be  found  in  Xu’s  thesis,  [17]. 

Also  a  paper  on  the  level  set  implementation  of  DSD  was  completed 
during  this  period  in  the  grant  and  it  was  later  published  in  the  Journal 
of  Computational  Physics,  [2].  The  level  set  code  used  for  this  paper  was 
DKAPPA2D,  written  by  Tariq  Aslam  as  part  of  his  PhD.  thesis,  [16]. 
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1.2.3  Summer  of  1995:  Installation  of  the  DSD  Wave  Tracker 
into  EPIC  Begins 

In  the  summer  of  1995,  Stewart,  Aslam  and  Directorate  engineer  Michael 
Nixon  started  the  process  of  installing  the  DSD  -  wave  tracker,  DKAPPA2D 
into  the  Eglin-based  code  EPIC,  A  similar  installation  for  MESA  had  already 
been  carried  out  in  collaboration  with  John  Bdzil  and  colleagues  at  LANL, 
again  with  the  use  of  DKAPPA2D.  In  addition  a  simple  DSD  program  burn 
algorithm  that  had  been  devised  by  Stewart  and  Bdzil  was  implemented  for 
the  polytropic  product  EOS.  The  June  of  1997  release  of  EPIC97  contained 
the  DSD  wavetracker  in  the  form  installed  by  Nixon. 

1.2.4  Summer  of  1995:  Beginnings  of  NONIDK  and  DSD  Char¬ 
acterization  of  Non-ideal  Explosives 

In  the  summer  of  1995,  Stewart  worked  out  the  basic  formulation  of  the 
simplest  DSD  theory  for  a  non-ideal  explosive  and  sketched  out  a  rough 
numerical  treatment  to  solve  the  nonlinear  eigenvalue  problem  for  the  — « 
response.  Later  this  effort  resulted  in  the  writing  of  NONIDK  a  code  that 
generates  the  Dn  —  n  response  curve  given  a  general  non-ideal  EOS  and 
reaction  rate  law  model.  This  early  work  was  documented  in  [9] . 

1.2.5  Summer  of  1995:  Stewart  and  Ruderman  Formulate  an 
Approach  to  Modeling  the  Mechanical  Ignition  of  Energetic 
Materials 

In  the  summer  of  1995,  while  visiting  Eglin  AFB,  Stewart  and  graduate  stu¬ 
dent  Ruderman  reviewed  some  of  the  approaches  being  used  by  Directorate 
engineers  to  understand  and  analyze  Taylor  Anvil  impact  experiments,  be¬ 
ing  carried  out  at  the  AWEFT  to  characterize  explosive  damage  and  igni¬ 
tion.  A  conclusion  was  drawn  that  the  state  of  modeling  justified  the  effort 
to  construct  a  new  first  principles  continuum  model.  The  explosive  HMX 
was  picked  as  a  target  material  for  the  new  modeling  effort  since  thermo¬ 
mechanical  data  for  HMX  is  relatively  abundant.  However,  the  modeling 
approach  to  be  devised  was  to  be  general.  It  was  decided  that  in  order  to 
properly  described  the  ignition  events  in  the  AWEFT  shear  and  impact  ex¬ 
periments,  that  the  phase  transformations  from  solid  to  liquid  and  liquid  to 
gas  and  subsequent  gaseous  combustion  all  needed  to  be  included.  If  pos¬ 
sible  the  physics  of  these  disparate  processes  should  be  combined  into  one 
self-consistent  continuum  model. 
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Ruderman  started  collecting  information  on  the  kinetics  and  properties 
and  phase  transformations,  particularly  melting  behavior  of  HMX,  Later 
that  fall  of  95  and  in  the  spring  of  1996  Stewart  and  Ruderman  began  a 
collaboration  with  a  UIUC  TAM  Professor  Eliot  Pried,  to  construct  a  phase- 
field  continuum  modeling  of  energetic  materials.  The  model  was  designed 
to  be  thermodynamically  consistent  in  three  dimensions.  The  model  has  an 
order  parameter  as  a  field  variable  that  is  used  to  describe  the  material  as 
a  solid,  liquid  and  gas.  Also  an  independent  reaction  progress  variable  was 
introduced. 

1.2.6  Fall  1995  through  Spring  1996:  Validation  of  DSD  through 
Direct  Comparison  with  DNS  via  AXS 

In  his  thesis  Aslam  [16]  demonstrated  conclusively  that  one  could  do  numer¬ 
ical  experiments  on  detonation  fiows  for  a  number  of  different  experiments  in 
different  geometries,  (converging  channel  wedge,  expanding  channel,  and  cir¬ 
cular  arc)  and  make  direct  measurement  of  the  detonation  dynamics.  Prom 
the  experimental  records  of  the  diffraction  and  convergence  experiments  it 
was  shown  that  one  could  successfully  construct  consistent  correlations  that 
describe  such  independent  experiments.  Prom  measurement  of  the  shock 
dynamics  observed  in  numerical  experiments  one  was  able  construct  a  rela¬ 
tionship  between  the  shock  acceleration,  the  normal  shock  velocity  Dn^ 
and  the  shock  curvature,  k,  relation  without  recourse  to  any  theory.  One 
could  then  independently  solve  the  DSD-like  evolution  equations  obtained 
from  these  measurements  and  compare  the  shock  locations  of  the  DSD-like 
prediction  and  the  DNS. 

The  comparisons  were  found  to  be  excellent  across  a  huge  range  of  det¬ 
onation  velocities,  curvature  and  shock  acceleration.  The  quality  of  the 
predictions  and  the  fact  that  this  very  different  methodology,  much  more 
akin  to  pure  experimental  work,  achieved  such  good  results  was  startling. 
This  work  also  indicated  the  extreme  importance  of  retaining  the  shock  ac¬ 
celeration  term  Dn  in  both  theory  and  computation  as  its  inclusion  made 
dramatic  improvements  in  the  quality  of  the  DSD  wave  front  predictions.  In 
1999  the  paper,  T.  Aslam  and  D.  S.  Stewart,  "Detonation  shock  dynamics 
and  comparison  with  direct  numerical  simulation” ,  appeared  in  Combustion 
Theory  and  Modeling,  [12]. 
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1.2.7  Summer  of  1996:  The  Second  AFOSR  Grant,  F49620-96- 
1-0260  Starts 

In  June  of  1996,  a  second  AFOSR  Grant,  F49620-96-1-0260  funded  by  Dr. 
Arje  Nachman,  or  AFOSR  Physical  Mathematics  Directorate,  entitled  ’’An¬ 
alytical  Investigations  on  Detonation  Theory  and  Mechanical  Ignition  of 
Condensed  Explosives”  started.  The  topical  coverage  of  the  grant  was  de¬ 
signed  to  augment  and  complement  the  research  goals  of  this  grant,  funded 
by  the  Directorate.  The  AFOSR  grant  provided  funds  for  the  Salary  of  Dr. 
Mark  Short  as  a  Research  Associate.  Dr.  Short  published  a  number  of  pa¬ 
pers  with  Stewart  on  the  stability  of  detonation,  [4],  [10],  [13].  Computing 
resources  and  travel  expenses  were  made  available  from  this  grant  to  sup¬ 
port  Dr.  Short.  Dr.  Short  visited  Eglin  AFB  in  the  summer  of  1997  and 
1998  and  began  to  interact  with  engineers  at  the  AWEFT  facility  on  issue 
associated  with  ignition  of  energetic  materials. 

1.2.8  Summer  of  1997:  Formulation  of  DSD  -  Cut,  Explosive 
System’s  Detonation  Shock  Pressure  Maps 

A  visit  to  Eglin  AFB,  in  the  Summer  of  1997  by  Stewart  witnessed  a  num¬ 
ber  of  important  developments  for  enhanced  DSD-design.  In  particular,  in 
response  to  Directorate  research  on  adaptable  warheads,  there  was  a  discus¬ 
sion  about  using  DSD  in  new  ways  and  taking  more  advantage  of  the  codes 
previously  written.  Stewart  proposed  a  new  engineering  methodology,  sub¬ 
sequently  dubbed  DSD-cut  to  distinguish  it  from  the  different  use  of  DSD 
that  was  implemented  in  EPIC  and  MESA,  subsequently  named  DSD-push. 

The  object  of  DSD-cut  was  to  use  DSD-computed  estimates  of  the  reac¬ 
tion  zone  structure  of  nonideal  explosives  to  compute  the  peak  detonation 
shock  pressures  and  the  reaction  zone  profiles  generated  by  detonation  prior 
to  the  impact  on  a  metal  target.  Specifically  the  focus  was  on  cutting  met¬ 
als.  This  work  was  deemed  particularly  relevant  for  explosives  (typical  of 
explosives  used  in  insensitive  munitions)  with  reaction  zones  that  may  be  a 
long  as  1/3  the  thickness  of  the  metal  target  plate. 

The  engineering  strategy  and  a  list  of  procedures  and  codes  needed  to 
carry  out  the  DSD-engineering  design  calculations  was  outlined  later  in  a 
TAM  technical  report  no  869,  ”  Detonation  Shock  Dynamics:  Application 
for  precision  cutting  of  metal  with  detonation  waves”,  [7].  Importantly  a 
very  simple  concept  of  mapping  the  detonation  shock  pressure  in  the  inte¬ 
rior  of  an  explosive  is  elucidated.  Shock  pressure  maps  allows  a  designer 
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immediate  feedback  on  the  placement  of  detonator  points  via  the  consistent 
use  of  the  DKAPPA2D  and  DKAPPA3D  class  of  codes.  NOIDK  provided  a 
means  to  generate  a  new  Dn  —  k  relations  for  an  explosive  used  in  a  design. 

Figures  1.  2.  and  3.  show  an  expample  of  the  DSD-based  detonation 
shock  pressure  maps  computed  in  the  interior  of  an  simple  shaped  charge 
geometry,  taken  from  [7]. 

1.2.9  Summer  1997;  Ehancements  to  NONIDK  to  Include  Shock 
Acceleration  Effects 

During  the  summer  of  1997,  Stewart  and  Nixon  worked  to  make  NONIDK 
platform  independent.  It  was  necessary  to  replace  the  numerical  integration 
subroutine  in  NONIDK  so  that  it  could  run  on  Eglin  platforms.  More  impor¬ 
tantly  Stewart  added  shock  acceleration  effects  into  a  formulation  of  DSD 
theory  that  could  be  used  for  condensed  explosives  described  by  a  nonideal 
EOS,  This  work  led  to  a  wholly  new  DSD  -  characterization  code  dubbed 
NONIDK-hPLUS  (to  be  distinguished  from  the  earlier  NONIDK), 

However  at  about  the  same  time,  the  simple  (interpolated  gas-solid)  EOS 
treatments  used  to  represent  the  reactive  equation  of  state  (an  interpolated 
gas-solid  EOS  used  in  the  first  versions  of  NONIDK)  became  suspect.  It  was 
found  that  certain  Hugoniots  crossed  and  hence  a  new  effort  on  obtaining 
better  EOS  models  was  deemed  as  essential  in  order  to  accurately  model 
explosives, 

1.2.10  Fall  1997:  The  Start  of  the  Davis  EOS  Model 

In  November  of  1997,  W,  C,  Davis  visited  UIUC  specifically  to  help  with  en¬ 
hancing  the  reactive  EOS.  This  was  an  important  effort  since  the  continued 
work  on  DSD  -  characterization  of  non-ideal  explosives  depended  on  mak¬ 
ing  improvements  to  the  EOS  models.  The  advantage  of  the  Davis/UIUC 
EOS  is  that  the  EOS  forms  used  are  based  on  explicit  functional  forms  that 
could  be  used  directly  in  NONIDK  and  NONIDK-fPLUS.  The  most  recent 
version  of  this  combined  EOS  improvements  and  shock  acceleration  effects 
was  documented  in  [14]. 

1.2.11  Spring  1998:  DSD  Literature  Search,  Invited  Topical  Re¬ 
view  on  DSD  written 

In  the  Spring  of  1998  a  substantial  effort  was  made  by  Stewart  to  carry  out 
a  comprehensive  literature  review  on  shock  dynamics  and  in  preparation  of 
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Figure  1:  Design  sketch  for  the  worked  example  in  the  report  [7].  The 
explosive  is  PBX9502,  the  side  liners  are  steel  and  the  top  cup  liner  is  copper. 
The  dimensions  are  marked  as  shown. 
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Figure  2:  Grey-scale  and  line  contour  plot  of  the  detonation  shock  pressure 
map  for  one  detonator  placement  from  [7], 


book  on  detonation  theory.  An  invited  topical  review  was  written  for  the 
Proceedings  of  the  Combustion  Institute,  entitled,  ”  The  shock  dynamics  of 
multi-dimensional  condensed  and  gas  phase  detonations”,  [11], 

1.2.12  Spring  1998  -  Summer  1998;  Design  of  GIBBSID  by  Ru- 
derman  and  Stewart 

In  the  early  Spring  of  1998  Stewart  and  Ruderman  outlined  the  design  of  a 
new  code  specifically  designed  to  solve  the  uni-directional  motions  of  RSF- 
model.  This  code  infrastructure  was  planned  to  solve  one-dimensional,  time- 
dependent  PDEs  that  could  handle  the  multi-scale  physics  of  the  phase 
transitions  and  the  ignition  events  in  the  HMX  model  The  code  suite  was 
dubbed  GIBBSID. 

It  was  also  decided  to  use  a  simplified  treatment  of  the  material  bound¬ 
ary  where  the  displacement  could  be  input.  It  was  also  decided  to  use 
ENO  schemes  and  to  adopt  some  of  the  same  numerical  strategies  found 
in  AXS.  Ruderman  began  the  code  development  with  extensive  testing  and 
code  validation.  A  copy  of  Ruderman’s  computer  program  used  to  simula¬ 
tion  phase  transition,  (called  GIBBSID)  was  forwarded  to  Eglin  engineers 
Michael  Nixon  and  K.  Vanden.  This  work  was  recorded  in  ’’Thermomechan¬ 
ical  Modeling  of  Energetic  Materials”,  [15]  and  was  presented  at  the  10th 
Detonation  Symposium  in  Snowmass  CO  in  August  1998.  A  more  detailed 
version  of  the  work  can  be  found  in  Ruderman’s  thesis,  [19].  A  substan¬ 
tial  revision  of  this  work  is  was  being  prepared  for  publication  with  current 
UIUC  graduate  student  Jack  Yoh. 

1.2.13  Fall  1998:  Preliminary  DSD-Characterization  of  PBXN9 
with  NONIDK-fPLUS 

In  the  Fall  of  1998,  the  Davis  EOS  was  implemented  in  NONIDK 4-PLUS. 
A  preliminary  characterization  of  PBXN9  was  carried  out  and  rate  law  was 
modified  so  that  the  DSD,  Dn  —  «  curve  for  zero  shock  acceleration  matched 
experimental  data  of  L.  Hull  of  LANL.  Estimates  of  the  shock  acceleration 
behavior  (contours  of  constant  Dn)  ior  PBXN9  were  generated  automatically 
from  NONIDK-f-PLUS.  This  preliminary  DSD  characterization  of  PBXN9 
was  reported  to  Directorate  personnel  in  a  March  1999  briefing  at  UIUC. 
The  work  is  also  recorded  in  [14]. 
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1.2.14  Winter  -  Spring  1999:  Simplified  Derivation  of  the  Theory 
of  DSD 

While  working  on  draft  introductory  chapters  of  the  book,  Stewart  sought 
to  simplify  the  presentation  of  the  DSD  theory.  In  one  of  the  drafts  Stewart 
used  the  new  formulation  originally  developed  in  the  topical  review  paper, 

[11]  to  re-calculate  the  ~  relation  first  given  in  [8],  With  Yaos 

help,  the  new  presentation  is  simpler  and  contains  all  the  previous  results. 
More  importantly  it  generates  new  and  simpler  analytic  formulas  for  the 
ideal  EOS  case.  The  ideal  EOS  case  is  an  important  test  case  and  has  been 
used  to  validate  NONIDK-l-PLUS.  The  new  asymptotics  are  also  compared 
favorably  against  the  numerical  experiments  of  Aslam  and  Stewart  found  in 

[12] ,  One  now  has  asymptotics,  a  general  numerical  treatment  and  direct 
numerical  simulation  all  in  agreement. 

1.2.15  Spring  1999:  Nineteen  Point  Detonation  Array  Simula¬ 
tions 

In  April  1999,  Stewart  visited  Directorate  personnel  at  Eglin  AFB  at  AWEFT. 
After  consultation  with  Dr.  Joseph  Foster,  and  Dr.  David  Lambert,  a  3- 
dimensional  DSD  wave  tracker  simulation  was  run  at  UIUC  for  a  19  pt 
hexagonal  pattern  array  of  detonators,  for  the  explosive  PBXN9  simulant, 
inside  a  metal  cylinder  The  tools  that  were  used  are  DKAPPA3D  and  NON- 
IDK-hPLUS  and  the  design  logic  discussed  in  [7],  The  Dn  =  0  curve  is  used 
as  the  Dn  —  k  relation  for  the  explosive.  A  version  of  DKAPPA3D  that  was 
set-up  to  carry  out  the  19  point  detonator  array  was  transferred  to  Eglin 
engineer  Michael  Nixon,  Figure  4.  shows  breakout  time  one  millimeter  from 
the  top  surface  of  a  19  point  detonator  system,  A  corresponding  detona¬ 
tion  shock  pressure  map  of  the  top  surface  at  break  out  is  easily  computed 
(not  shown).  The  computations  were  done  with  DKAPPA3D,  by  Stewart 
and  Yao  in  collaboration  with  Eglin,  AWEFT  egineers  David  Lambert  and 
Joseph  Foster.  The  Dn  —  n  relation  used  was  for  a  simulant  of  PBXN9. 


2  Archival  papers,  Ph.D.  theses  and  Reports 

Below  we  list  the  most  important  papers,  reports  and  Ph.D.  theses  that  ben¬ 
efited  wholly  or  in  part  from  the  resources  of  this  grant.  The  entries  appear 
in  more  or  less  chronological  order  and  the  reference  number  corresponds  to 
the  text  of  the  preceding  section.  An  abstract  or  summary  of  each  entry  is 
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Figure  4:  Map  of  the  arrival  time  of  a  19  point  detonator  system  computed 
with  DKAPPA3D  for  a  PBXN9  explosive  simulant,  May  1999.  The  arrival 
time  map  corresponds  to  the  cut  pattern  of  a  top  metal  plate  induced  by 
detonation  in  the  explosive  initiated  by  the  19  point  array. 
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included  to  summarize  the  research  findings  in  more  detail. 
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intended  application,  as  opposed  to  a  finite  volume  implementation.  Our 
solver  allows  for  the  placement  of  internal  reflective  boundaries  and  the 
standard  inflow  and  outflow  and  reflective  boundaries  at  the  edge  of  the 
domain.  We  discuss  the  spatial  discretization  and  the  temporal  integration 
schemes,  upwinding  and  flux  splitting  and  the  combined  use  of  the  Lax- 
Friedrichs  and  Roe  schemes  to  solve  for  the  required  fluxes.  A  complete 
description  of  the  pointwise  internal  boundary  method  is  given.  An  overall 
summary  of  a  representative  code  structure  is  given.  We  provide  details 
on  the  verification  of  our  integrated  set  of  algorithms  that  resulted  in  an 
application  code.  We  demonstrate  the  order  of  convergence  for  test  prob¬ 
lems.  Two  example  applications  from  measurement  of  detonation  shock 
dynamics  and  deflagration  to  detonation  transition  in  porous  energetic 
materials  are  presented. 

[7]  Stewart,  D.  S,  Detonation  shock  dynamics:  Application  for  precision 
cutting  of  metal  with  detonation  waves,  TAM  Report,  No.  869,  UILU- 
ENG-97-6030,  October  1997. 

Abstract.  This  technical  report  starts  with  an  introduction  about  the 
engineering  use  of  detonation  shock  dynamics  and  its  advantages.  Next 
there  is  a  review  of  the  relevant  basics  of  DSD-theory  that  include  the  re¬ 
duced  governing  equations,  the  Rankine-Hugoniot  relations  and  CJ-state, 
the  formulation  of  the  basic  eigenvalue  problem  that  is  used  to  relate  the 
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normal  detonation  velocity  and  the  shock  curvature  and  representative 
results  for  PBX9502.  The  next  section  explains  how  to  make  shock  pres¬ 
sure  maps  in  the  interior  of  an  explosive  system.  First  the  DSD  wave 
tracker  is  run  on  the  explosive  domain,  then  the  shock  pressure  is  com¬ 
puted,  Comments  are  made  on  the  interpretation  of  the  maps  for  design 
purposes. 

[8]  Buckmaster,  J.,  Short  M.  and  Stewart,  D.  S.  The  use  of  activation  energy 
asymptotics  in  detonation  theory,  with  Comment  on  ’’Multidimensional 
stability  analysis  of  overdriven  gaseous  detonation,  Phys.  Fluids  9,  3764, 
Physics  of  Fluids,  10,  3027-3030  (1998). 

Abstract.  A  powerful  defense  is  given  for  the  use  of  activation  energy 
asymptotics  in  opposition  to  erroneous  and  misleading  work  carried  out 
by  Clavin  and  Williams. 

[9]  Stewart,  D.  S.  and  Yao,  J.  The  normal  shock  velocity-curvature  relation¬ 
ship  for  materials  with  non-ideal  equation  of  state  and  multiple  turning 
points.  Combustion,  113,  224-235  (1998). 

Abstract.  We  present  a  model  and  simple  to  implement  numerical  pro¬ 
cedure  that  obtains  the  normal  detonation  shock  velocity  curvature  rela¬ 
tionships  for  an  explosive  material  with  non-ideal  equation  of  state  and 
an  arbitrary  reaction  rate  law.  In  addition  we  illustrate  numerically  (for 
a  non-ideal  equation  of  state)  and  analytically  (for  an  ideal  equation  of 
state  with  a  large  activation  energy  rate  law)  that  for  sufficient  rate-state 
-sensitive  explosives,  the  response  curve  can  have  two  turning  points  such 
that  the  curve  has  a  Z-shape.  The  top  branch  of  the  Z  response  curve 
has  been  previously  been  associated  with  detonation  extinction  at  a  crit¬ 
ical  curvature.  The  bottom  branch  can  be  possibly  associated  with  low 
velocity  detonation  and  rapid  change  from  low  order  detonation  to  high 
order, 

[10]  Short,  M.  and  Stewart,  D.S.  Cellular  detonation  stability:  A  normal 
mode  linear  analysis.  Journal  of  Fluid  Mechanics,  368,  229-262  (1998). 

Abstract.  A  detailed  investigation  of  the  hydrodynamic  stability  of  a 
steady,  one-dimensional  detonation  to  transverse  linear  disturbances  in  an 
ideal  gas  undergoing  an  irreversible,  unimolecular  reaction  with  an  Arrhe¬ 
nius  rate  constant  is  conducted  via  a  normal-mode  analysis.  The  method 
of  solution  is  an  iterative  shooting  method,  which  integrates  between  the 
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detonation  shock  and  reaction  equilibrium  point.  Variations  in  the  distur¬ 
bance  growth  rates  and  frequencies  with  transverse  wavenumber,  together 
with  two-dimensional  neutral  stability  curves  and  boundaries  for  all  un¬ 
stable  low-  and  high-frequency  modes  are  obtained  for  varying  detonation 
bifurcation  parameters.  These  include  the  detonation  overdrive,  chemical 
heat  release  and  reaction  activation  energy.  Spatial  perturbation  eigen¬ 
function  behavior  and  phase  and  group  velocities  are  also  obtained  for 
selected  sets  of  unstable  modes.  Results  are  presented  for  both  Chapman- 
Jouguet  and  overdriven  detonation  velocities.  Comparisons  between  the 
earlier  pointwise  determination  of  stability  and  interpolated  neutral  stabil¬ 
ity  boundaries  obtained  by  Erpenbeck  are  made.  Possible  physical  mecha¬ 
nisms  which  govern  the  wavenumber  selection  underlying  the  initial  onset 
of  either  regular  or  irregular  cell  patterns  are  also  discussed. 

[11]  Stewart,  D.  S.  The  shock  dynamics  of  multi-dimensional  condensed  and 
gas  phase  detonations.  Proceedings  of  the  27th  International  Symposium 
on  Combustion,  pp.  2189-2205  (1998) 

Abstract,  Detonations  are  comprised  of  broad  detonation  shock  sup¬ 
ported  by  thin  reaction  zone.  Approximation  base  on  weak  shock  curva¬ 
ture  measured  on  the  inverse  reaction  zone  scale,  and  quasi-steady  flow, 
measured  on  the  particle  passage  time  through  the  reaction  zone  can 
be  used  to  simplify  the  mathematical  description  of  detonations  that  are 
governed  by  the  gasdynamic  equations  for  a  reacting  flow.  When  the  deto¬ 
nation  reaction  zone  contains  a  sonic  locus  it  is  possible  to  derive  intrinsic 
(coordinate  independent)  partial  differential  equation  for  the  lead  deto¬ 
nation  shock’s  motion  in  terns  of  the  normal  detonation  shock  velocity, 
the  shock  curvature  and  higher  normal  time  derivatives.  We  refer  to  this 
collection  of  theory  and  supporting  experimental  results  as  Detonation 
Shock  Dynamics  after  Whitham’s  Geometrical  Shock  Dynamics.  The  re¬ 
duce  detonation  dynamics  is  based  on  the  concept  of  an  eigenvalue  (sonic) 
detonation,  an  idea  that  goes  back  to  the  original  investigations  in  the 
1940’s.  We  present  a  review  of  the  theoretical  and  experimental  develop¬ 
ments  and  attempt  to  update  Pickett  and  Davis’  discussion  of  work  prior 
to  1980.  We  give  examples  of  the  theory  and  applications  which  include: 
i)  weakly-curved,  quasi-steady,  near-CJ  detonation  ii)  critical  detonation 
curvature  iii)  quasi-steady  extinction  and  ignition  (and  low  velocity  det¬ 
onation)  iv)  shock  acceleration  effects  and  v)  cellular  and  pulsating  det¬ 
onation  in  gases.  We  also  review  the  engineering  method  of  Detonation 
Shock  Dynamics  as  it  is  applied  to  explosive  systems. 
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[12]  Aslam,  T.D.  and  Stewart,  D.S.  Detonation  shock  dynamics  and  com¬ 
parisons  with  direct  numerical  simulation,  Combustion  Theory  and  Mod¬ 
eling,  3,  77-101  (1999) 

Abstract.  Comparisons  between  direct  numerical  simulation  (DNS)  of 
detonation  and  detonation  shock  dynamics  (DSD)  is  made.  The  theory 
of  DSD  defines  the  motion  of  the  detonation  shock  in  terms  of  the  in¬ 
trinsic  geometry  of  the  shock  surface,  in  particular  for  condensed  phase 
explosives  the  shock  nomrl  veloticy,  the  normal  acceleration,  and 
the  total  curavture,  k.  In  particular,  the  properties  of  three  intrinsic  front 
evolution  laws  are  studied  and  compared.  This  are  (i)  constant  speed  det¬ 
onation  (Huygens’s  construction),  (ii)  curvature-dependent  speed  propa¬ 
gation  {Dn  —  K  -  relation)  and  (iii)  curvature  speed-dependent  acceleration 
{Dn  —  Dn  —  K  relation).  We  show  that  it  is  possible  to  measure  shock  dy¬ 
namics  directly  from  simulation  of  the  reactive  Euler  equations  and  that 
subsequent  numerical  solution  of  the  intrinsic  partial  differetial  equation 
for  the  shock  motion  (e.g.  a  Dn—Dn—K  relation)  reproduces  the  computed 
shock  motion  with  high  precision. 

[13]  Short,  M.  and  Stewart,  D.  S.  The  multi-dimensional  stability  of  weak 
heat  release  detonations.  J.  Fluid  Mech.,  s382,  109-135  (1999). 

Abstract.  The  stability  of  an  overdriven  planar  detonation  wave  is  exam¬ 
ined  for  a  one-step  Arrhenius  reaction  model  with  an  order  one  post-shock 
temperature-scaled  activation  energy  in  the  limit  of  a  small  post-shock 
temperature-scaled  heat  release.  The  ratio  of  specific  heats,  is  taken  to  be 
close  to  one.  Under  these  assumptions,  which  cover  a  wide  range  of  realis¬ 
tic  physical  situations,  the  steady  detonation  structure  can  be  evaluated 
explicitly,  with  the  reactant  mass  fraction  described  by  an  exponentially 
decaying  function.  The  analytical  representation  of  the  steady  structure 
allows  a  normal-mode  description  of  the  stability  behavior  to  be  obtained 
via  a  two-term  asymptotic  expansion  in  the  heat-release.  The  resulting 
dispersion  relation  predicts  that  for  a  finite  overdrive  the  detonation  is 
always  stable  to  two-dimensional  disturbances.  For  large  overdrives,  the 
identification  of  regimes  of  stability  or  instability  is  found  to  depend  on  a 
choice  of  distinguished  limit  between  the  heat  release  and  the  detonation 
propagation  Mach  number.  Regimes  of  instability  are  found  to  be  charac¬ 
terized  by  the  presence  of  a  single  unstable  oscillatory  mode  over  a  finite 
range  of  wavenumbers. 
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[14]  Stewart,  D.  S.  and  Yao,  J.,  Computation  of  shock  acceleration  effects  on 
detonation  shock  dynamics  for  explosives  described  by  general  equation 
of  state,  to  appear  Proceedings  of  the  28th  International  Symposium  on 
Combustion. 

Abstract.  The  inclusion  of  detonation  shock  acceleration  effects  leads 
to  an  extended  theory  of  Detonation  Shock  Dynamics  (DSD).  The  shock 
motion  is  described  by  an  intrinsic  partial  differential  equation  specified 
in  terms  of  the  normal  shock  velocity,  the  normal  shock  acceleration 
Dri  and  the  curvature  /c.  Earlier  developments  were  based  on  analytical 
(asymptotic)  calculations  and  carried  out  for  the  polytropic  equation  of 
state  that  made  detailed  analysis  tractable.  But  the  demands  of  quanti¬ 
tative  accuracy  for  engineering  design  require  that  real  equation  of  state 
constitutive  forms  be  used  with  more  functional  complexity.  In  this  pa¬ 
per  we  present  a  numerical  approach  that  can  be  used  to  compute  the 
Dn,  -Dn,  /^-relation  for  general  EOS  and  rate  law  forms. 

[15]  Ruderman,  G.,  Stewart  D.  S.  and  Fried,  E.  Modeling  the  Mechani¬ 
cal  Ignition  Energetic  Materials,  Proceedings  of  the  11th  (International) 
Symposium  on  Detonation,  Office  of  Naval  Research,  pp  554-561  (2000) 

Abstract,  A  new  theory  of  continuum  thermomechanics  involving  mi¬ 
croforces,  or  forces  which  drive  the  evolution  of  the  microstructure  in  a 
material,  is  applied  to  the  modeling  of  energetic  materials,  specifically, 
HMX.  Using  this  theory,  a  set  of  partial  differential  equation  are  de¬ 
rived  in  a  consistent  thermodynamic  framework  which  models  HMX  as 
a  nonlinearly  viscoelastic  solid  capable  of  undergoing  phase  transition  to 
a  viscous  liquid  and  ideal  gas,  as  well  as  a  combustion  process  from  un¬ 
reacted  HMX  to  reacted  products.  This  system  was  then  simplified  to  a 
set  of  one-dimensional  model  problems,  which  were  solved  numerically  us¬ 
ing  advance  computation  techniques  including  essentially  non-oscillatory 
(ENO)  methods. 

2.2  Phd.  Theses  Supported  by  the  Grant 

[16]  Aslam,  T.  D.,  Investigations  of  Detonation  Shock  Dynamics,  Ph.  D. 
thesis,  University  of  Illinios,  1996,  Theoretical  and  Applied  Mechanics, 
thesis  advisor  D.  S.  Stewart 

Abstract.  A  detonation  is  a  combustion  driven  shock  wave.  Typically, 
a  detonation  will  consist  of  an  inert  shock  followed  by  a  region  of  chem- 
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ical  reaction  referred  to  as  the  reaction  zone.  Detonations  have  a  wide 
variety  of  engineering  applications,  from  obvious  military  use  to  explosive 
welding,  hard  rock  mining  and  materials  processing.  Detonations  can  oc¬ 
cur  in  a  variety  of  materials,  including  gases  (such  as  premixed  hydrogen 
and  oxygen),  liquids,  and  solid  explosives.  Of  particular  interest  in  det¬ 
onation  problems  is  the  motion  of  the  detonation  shock.  Changes  to  the 
reaction  zone  may  cause  large  variations  in  the  strength  and  speed  of  the 
detonation  front,  so  it  can  not  be  ignored  in  modeling  detonations.  For 
typical  explosives,  the  reaction  zone  may  be  thousands  of  times  smaller 
than  the  engineering  scale.  This  multi-scaled  nature  of  detonation  can 
pose  problems  when  trying  to  predict  the  motion  of  the  detonation  front. 

Detonation  shock  dynamics  is  an  asymptotic  theory  whose  key  result  is  an 
intrinsic  partial  differential  equation  for  the  dynamics  of  the  detonation 
shock  front.  It  will  be  demonstrated  that  the  theory  can  predict  several  as¬ 
pects  of  unsteady  multi-dimensional  detonation  accurately.  Three  intrinsic 
relations  will  be  examined  and  compared  with  direct  numerical  simula¬ 
tion.  Their  relevance  to  modeling  detnation  dynamics  will  also  be  given. 
Numrical  methods  based  on  level-set  ideas,  will  be  given  for  propagating 
multi-dimensional  detonation  fronts  in  arbitrarily  complex  geometries. 

[17]  Xu,  Shaojie,  Modeling  and  numerical  simulation  of  deflagration  to  det¬ 
onation  transition  in  porous  energetic  materials,  Ph.  D.  thesis,  University 
of  Illinios,  1996,  Theoretical  and  Applied  Mechanics,  thesis  advisor  D.  S. 
Stewart 

Abstract.  An  understanding  of  the  deflagration- to-detonation  transition 
(DDT)  in  porous  energetic  materials  is  important  to  various  engineering 
applications.  Safety  issues  of  explosives,  specifically  the  prevention  of  ac¬ 
cidental  detonation  is  damaged  explosive,  is  one  example.  The  transition 
process  is  complex  in  nature  and  involves  both  subsonic  and  supersonic 
phenomena  such  as  compaction  and  detonation  respectively.  The  presence 
of  multiple  time  and  length  scales  makes  the  study  of  DDT  challenging  and 
interesting.  Modern  computational  technology  has  made  a  great  impact 
on  the  study.  Numerical  simulation  becomes  an  important  a  necessary 
means  to  carry  out  the  study,  especially  for  problems  arising  from  engi¬ 
neering  applications.  Two  topics  related  to  multi-dimensional  simulation 
of  DDT  in  energetic  materials  are  presented  in  this  work. 

The  objective  of  the  first  part  is  to  develop  a  model  that  has  a  relatively 
simple  mathematical  structure.  Two-phase  mixture  theory  has  been  used 
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over  the  past  two  decades  to  describe  DDT  in  porous  energetic  materials. 
Conservation  laws  of  mass,  momentum  and  energy  are  written  indepen¬ 
dently  for  two  materials,  and  the  interactions  between  the  phases  are  rep¬ 
resented  by  source  terms.  Models  constructed  by  the  theory  usually  have 
complicated  structures  represented  by  a  large  set  of  governing  equations 
and  complex  wave  structures.  Hence,  it  is  difficult  to  resolve  them  nu¬ 
merically  with  the  computational  power  available.  Three  simplified  DDT 
models  name  Bdzil-Kapila-Stewart  (BKS),  solid- void-gas  (SVG)  and  gas- 
interpolated  solid  Stewart  Prasad  Asay  (GISPA),  are  considered  in  the 
current  study;  among  them  SVG  and  GISPA  are  newly  developed.  These 
models  are  base  on  a  single  velocity  formulation  and  thus  have  fewer 
governing  equations  and  simple  wave  structures.  A  comparative  study  is 
carried  out  based  on  the  three  models.  The  purpose  of  the  study  is  to 
identify  a  minimum  set  of  elements  possessed  by  a  model  to  make  predic¬ 
tions  consistent  with  experiments.  The  study  shows  that  the  SVG  mode  is 
the  best  among  the  three  models  over  a  wide  range  of  experimental  data. 

In  addition  to  model  simplification,  effort  has  been  focused  on  chemical 
reaction  kinetic  modeling,  which  describes  energy  release  during  the  re¬ 
action  process.  Previously  used  chemical  kinetic  models,  or  reaction-rate 
laws,  are  not  adequate  to  predict  detonation  initiation.  A  new  reaction 
rate  law  is  developed  which  employs  two  different  rate  functions  associ¬ 
ated  with  slow  and  fast  energy-release  processes  that  have  been  observed 
in  experiments.  With  the  new  rate  law,  the  SVG  and  GISPA  models  are 
able  to  predict  qualitatively  and  quantitatively  all  the  events  recorded  in 
one-dimensional  DDT  tube  experiments. 

The  second  part  of  the  study  is  devoted  to  the  development  of  a  high- 
quality  numerical  method,  suitable  for  multi-dimensional  direct  numeri¬ 
cal  simulation  of  DDT  with  non  -linear  material  on  a  complex  geometry. 
The  method  is  developed  through  an  integration  of  contemporary  shock¬ 
capturing  methods.  The  new  finite-difference  scheme  adopts  the  method- 
of-lines  approach,  which  allows  for  independent  temporal  and  spatial  dis¬ 
cretizations.  The  temporal  integration  uses  a  third-order  Runge-Kutta 
method  with  the  property  of  total  variation  diminishing  (TVD),  while 
the  spatial  integration  employs  a  fourth-order  essentially  non-oscillatory 
(ENO)  scheme.  In  order  to  implement  the  scheme  for  nonlinear  mate¬ 
rial,  a  Roe’s  linearized  Riemann  solver  is  developed  in  two  dimension 
for  an  equation  of  state  that  describe  an  HMX  material  To  fulfill  the 
needs  of  engineering  applications,  an  internal  boundary  method  is  de- 
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veloped  on  a  structured  grid.  The  internal  boundary  algorithm  allows  a 
two-dimensional  non-deformable  body  of  arbitrary  shape  to  be  inserted 
in  a  flow  field.  A  second-order  reflective  boundary  condition  is  implement 
along  the  internal  boundary.  The  new  scheme,  equipped  with  the  internal 
boundary  algorithm,  provides  great  flexibility  for  numerical  simulation  of 
engineering  problems  with  complex  boundaries.  A  two-dimensional  DDT 
simulation  is  carried  out  to  simulate  a  blunt  body  impact  on  energetic 
materials.  A  bullet  with  a  circular  head  is  imbedded  in  a  square  body 
of  HMX  material  and  is  subjected  to  an  impulsive  motion.  The  radius 
effect  on  transition  to  detonation  is  studied.  A  sharp  body  impact  is  also 
considered  where  the  circular  bullet  head  is  replaced  by  a  wedge. 

[18]  Yao,  J.  The  dynamics  of  multi-dimensional  detonation.  Ph.  D.  thesis, 
University  of  Illinois,  1996,  Theoretical  and  Applied  Mechanics,  thesis 
advisor  D,  S.  Stewart 

Abstract.  An  asymptotic  theory  is  presented  for  the  dynamics  of  det¬ 
onation  when  the  radius  of  curvature  of  the  detonation  shock  is  large 
compared  with  one-dimensional  steady  Chapman-Jouguet  (CJ)  deto¬ 
nation  reaction-zone  thickness.  The  analysis  considers  additional  time- 
dependence  in  the  slowly-varying  reaction  zone  than  that  considered  in 
previous  works.  The  detonation  is  assumed  to  have  a  sonic  point  in  the 
reaction-zone  structure  behind  the  shock,  and  is  referred  to  as  an  eigen¬ 
value  detonation.  A  new  iterative  method  is  used  to  calculate  the  eigen¬ 
value  relation,  which  ultimately  is  expressed  as  an  intrinsic  partial  dif¬ 
ferential  equation  (PDE)  for  the  motion  of  the  shock  surface.  Two  cases 
are  considered  for  an  ideal  equation  of  state.  The  first  corresponds  to 
a  model  of  a  condensed  phase  explosive,  with  model  reaction-rate  sen¬ 
sitivity,  and  the  intrinsic  shock  surface  PDE  is  a  relation  between  the 
normal  detonation  shock  velocity,  the  first  normal  time  derivative  of 
the  normal  shock  velocity,  Z)^,  and  the  shock  curvature  k.  The  second 
case  corresponds  to  a  gaseous  explosive  mixture,  with  the  large  reaction- 
rate  sensitivity  of  Arrhenius  kinetics,  and  the  intrinsic  shock  surface  PDE 
is  a  relation  between  the  normal  detonation  shock  velocity  its  first 
and  second  normal  time  derivatives  I)^,  the  shock  curvature  and 
the  first  normal  time  derivative  of  the  curvature  k.  For  the  second  case, 
one  obtains  a  one-dimensional  theory  of  pulsation  of  plane  G  J  detonation 
and  a  theory  that  predicts  the  evolution  of  self-sustained  cellular  detona¬ 
tion.  Versions  of  the  theory  include  the  limit  of  near-CJ  denotation,  and 
the  limit  in  which  the  normal  detonation  velocity  is  significantly  below 
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its  CJ  value.  The  curvature  of  the  detonation  can  also  be  of  either  sign 
corresponding  to  either  diverging  or  converging  geometry. 

The  linear  instability  of  a  weakly  curved  slowly  varying  detonation  wave  is 
also  investigated  under  the  assumption  of  frozen  curvature.  The  governing 
equations  and  the  boundary  conditions  required  to  formulate  the  instabil¬ 
ity  problem  are  derived.  The  steady  Dn  —  k  relation  and  the  quasi-steady 
state  of  the  weakly  curved  detonation  have  been  obtained  numerically. 
The  eignvalues  of  the  acoustic  instability  are  calculated  by  a  numerical 
shooting  method. 

[19]  Ruderman,  G.  A.,  A  Continuum  Thermomechanical  Model  for  Ener¬ 
getic  Materials,  Ph.  D.  thesis,  University  of  Illinois,  1998,  thesis  advisor 
D.  S.  Stewart,  TAM. 

Abstract.  Thermomechanical  modeling  of  energetic  materials,  for  exam¬ 
ple  solid  rock  motor  propellants  and  explosives,  is  a  complex  problem  due 
to  the  large  number  of  behaviors  such  a  material  may  exhibit.  Experi¬ 
ments  have  shown  that  these  materials  are  nonlinearly  viscoelastic,  and 
may  also  experience  plastic  flow  permanent  deformation)  phases  changes 
(melting  and  vaporization  processes),  and  combustion.  In  addition,  these 
phenomena  are  often  strongly  coupled,  making  modeling  very  difflcult. 
Compounding  this  difficulty  further,  reliable  experimental  data  on  the 
properties  of  these  types  of  materials  are  quite  scarce. 

Applying  advanced  tools  of  continuum  thermomechanics,  we  have  devel¬ 
oped  a  fully  three-dimensional  framework,  which,  in  the  most  general 
form,  is  able  to  model  all  the  mentioned  behaviors  of  energetic  materials. 
The  concept  of  a  balance  of  micro- forces,  which  drives  changes  in  material 
microstructure,  is  employed  to  generate  thermomechanically  consistent 
equations  of  evolutions  for  combustion  and  phase  transition. 

The  model  is  the  simplified  to  a  set  of  three  model  problems:  the 
constant-volume  thermal  explosion,  one-dimensional  shear  loading,  and 
one-dimensional  longitudinal  loading.  These  model  problems  were  solved 
numerically  using  essentially  non-oscillatory  and  total  variation  dimin¬ 
ishing  methods.  The  solutions  reveal  extremely  rich  behavior,  including 
complex  wave  phenomena,  strain  localization  phenomena,  and  changes  of 
material  phase. 
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